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Abstract . 


A method for calculating ensemble averages from magnetic 
field data is described. A data set comprising approximately 
sixteen months of nearly continuous ISEE-3 magnetic field data 
is used in this study. Individual subintervals of this data, 
ranging from 15 hours to 15.6 days comprise the ensemble. The 
sole condition for including each subinterval in the averaqes is 
the degree to which it represents a weakly time-stationary 
process. Averages obtained by this method are appropriate for a 
turbulence description of the interplanetary medium. The 
ensemble average coi ilation length obtained from all subinter- 
vals is found to be 4.9 x 10 11 cm. The average value of the 
variances of the magnetic field components are in the approxi- 
mate ratio 8:9:10, where the third component is the local mean 
field direction. The correlation lengths and variances are 
found to have a systematic variation with subinterval duration, 
reflecting the important role of low-frequency fluctuations in 
the interplanetary medium. 
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Infer Qfluqt i o n 


Studies of solar wind fluctuations have been carried out by 
a number of investigators C Coleman, 1966; Jokipii and Coleman, 
1968; Sari and Ness, 1969; Belcher and Davis, 1971; Barnes, 

1979; Burlaga, 1985; among others!. Most previous studies have 
sought to characterize specific types of solar wind flows or to 
give examples of specific phenomenological features of inter- 
planetary fluctuations. Usually the number of intervals 
analyzed in a given study has not been large. Comparison of 
studies is hampered because the criteria for selecting data 
intervals have varied according to the objectives of the 
Investigation, as have the data analysis procedures. Only a 
small fraction of the available interplanetary data has been 
studied statistically. A considerably sharper perspective of 
the state of interplanetary magnetic fluctuations would be 
available if a systematic procedure for extracting convergent 
statistical estimates were applied to large numbers of field 
samples. In essence, one would like to have an appropriate 
algorithm for constructing statistically valid ensembles of 
solar wind data sets. The properties of such ensembles could 
then be used to investigate physical problems ranging from 
magnetohydrodynamic (MHD) turbulence to cosmic ray propagation 
and modulation. 

He describe here a method for accomplishing this task and 
apply it to solar wind magnetic field data acquired at 1 AU. 
Accurate determination of quantities relevant to the development 
of a turbulence description of the solar wind is a primary 
motivation in the present study; however, the bulk of the 
results do not rely on specific dynamical models of the inter- 
planetary medium. The methods used are an extension of those 
employed previously by Matthaeus and Goldstein C1982a,b! to 
Investigate the MHD turbulence properties of the solar wind. 


- 3 - 



The essence of our approach is that reproducible turbulence 
properties of the Magnetic field are average properties of a 
statistical ensemble , and that a large collection of appropri- 
ately chosen analyses of solar wind data Intervals constitute, 
in effect, rather precise knowledge of that ensemble. 

In the following sections, some background is presented and 
the method is presented in detail. Subsequently, the data base 
used is described, and the results presented. Finally, implica- 
tions of the results are discussed and directions for future 
work presented. 


Background . 

Statistical methods are useful in describing physical 
systems that exhibit ^reproducible behavior. Such methods have 
been an essential ingredient in turbulence theory CTaylor, 1938; 
Batchelor, 1970; Orszag, 19773 and experiments CGrant, Stewart 
and Moilliet, 1962; Hinze, 19753 for most of this century. In 
the statistical perspective, individual samples of the dynamical 
variables, which may be, for example, experimentally determined 
values of fluid velocity in a turbulent fluid, are taken to be 
particular values of random variables, or "realizations". A 
probability distribution function, specifying the relative 
likelihood of occurrence of the various realizations, is assumed 
to exist when it is not calculable from first principles. The 
use of probability methods for a deterministic system is most 
easily made plausible by appealing to the notion of an "ensem- 
ble", which may be thought of CVan Kampen, 19813 as the collec- 
tion of all possible samples of a system, or equivalently a 
collection of a very large number of similarly prepared systems. 
Expectation values calculated from the probability distribution 
are identical to averages over the ensemble. In effect, the 
probability distribution is just the tabulation of the frequency 
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of occurrence of specific values of the random variables in the 
ensemble. 

If the solar wind is viewed as a turbulent medium, it is 
natural to apply this type of method to its study. Unfor- 
tunately, the current level if knowledge (see, e. g.. Whang 
C19773) about the probability distribution is not sufficient to 
base these investigations on the probability function itself. 
Consequently, most work has concentrated on statements about 
empirically determined statistical quantities, such as correla- 
tion times and lengths, spectra and correlation functions 
CJokippi, 19733. If accurately determined, these quantities are 
mathematically related to the probability distribution, but 
their validity as estimates of ensemble averages is rarely 
considered in specific cases. In all cases, the scope of 
applicability of a statistical measurement in the solar wind 
must be judged in accordance with the def ensibility of the data 
selection criteria, the rationale for the method of analysis, 
and the stability of the results when data intervals and 
analysis parameters are varied. At the same time, physically 
useful statistical results must preserve information that may be 
both model -dependent and sensitive to analysis procedures. 

The nature of the interplanetary medium renders the task of 
living up to these ideals difficult, particularly when real- 
world limitations on the number, duration, sampling rates and 
instrumental accuracy of spacecraft experiments are taken into 
account. The solar wind contains structure ranging over a wide 
dynamic range of both spatial and temporal scales. At the 
largest scales and lowest frequencies, the solar wind exhibits 
organization, including systematic variations cf many physical 
properties with heliocentric distance as well as the global 
Parker spiral magnetic structure. In addition, large scale 
features such as high-speed streams and sector structure 
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CHolzer, 19793 , interaction regions , corotating and transient 
flow patterns CDryer et al., 1976; Burlaga, 19833 Influence and 
organize significant volumes of the Interplanetary plasma. 

There are also levels of organization at Intermediate scales 
that might be thought of as "coherent structures", such as 
shocks and MHD waves and discontinuities. Finally, present in 
all plasma variables, are broad-band stochastic fluctuations 
extending roughly from the reciprocal of the correlation time to 
beyond the ion-cyclotron frequency. These structural character- 
istics are important ingredients in the physics we study, but 
they also make the signals that must be selected and interpreted 
much more difficult to deal with than unstructured random 
processes, for example, the random telegraph signal CRice, 

19543. 

Researchers have responded to these complications with a 
variety of selection and analysis methods. Most often, selec- 
tion of data has proceeded with the goal of studying a particu- 
lar class of solar wind behavior. In many cases, the scientific 
conclusions depend to some degree on the method used. 

For example "Alfvfenic" periods are probably the most widely 
studied CColeman, 1968; Jokippi and Coleman, 1968; Belcher and 
Davis, 19713 because of thei\' frequency of occurrence and the 
likelihood that they represent MHD behavior, as well as because 
of their potential influence on solar and galactic cosmic rays. 
However, the conclusion that Alfvfenic periods consist entirely 
of "outward travelling waves" CBelcher and Davis, 1971; 
Dobrowolny, Mangeney and Veltri, 19803 has been questioned. 
Matthaeus and Goldstein C1982a3 found intervals which displayed 
outward travelling Alfvfenic correlations becween velocity and 
magnetic field vectors at frequencies corresponding to the 
inertial range of MHD turbulence, but correlations corresponding 
to "inwa •'■^-travelling" waves were found at lower frequencies. To 
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see this reversal of the sign of the cross hellcity, the data 
Interval must be sufficiently long. In addition, in one 
interval obtained near 5 AU, Alfv6nic fluctuations with periods 
within the inertial range were observed to propagate randomly 
both inward and outward. The characteristic mixed sign of the 
velocity-magnetic field correlation was also seen in two- 
dimensional MHD turbulence simulations CMatthaeus, Goldstein and 
Montgomery, 19833, which led to the suggestion that Alfvfenic 
periods may be produced by in situ MHD turbulence rather than 
representing a remnant of coronal processes. 

The processing of data can also influence statistical 
conclusions. The standard time-series analysis technique of 
time domain windowing COtncs and Enochsen, 19723 is justified by 
citing the requirement of strict periodicity when fast Fourier 
transform spectral analysis is used. However, windowing 
systematically reduces the complement of low frequency power, 
causing a reduction in estimates of correlation lengths and 
times CMatthaeus and Goldstein, 1982a3. Another type of data 
preprocessing which gives results that must be interpreted with 
care is exemplified by the analysis of Pioneer 6 data by Fisk 
and Sari C13733 (also see Sari and Valley C19763). They removed 
tangential discontinuities from selected data intervals before 
analysis with the goal of extracting properties of the "under- 
lying" random signal. This method is of questionable applicabi- 
lity to a turbulence description of the solar wind since MHD 
turbulence seems to have a propensity for producing nearly 
discontinuous vorticity and current structures CMatthaeus and 
Montgomery, 1980; Matthaeus, 19823. 

Data selection and processing can also Influence the 
sensitive issue of defining the mean magnetic field. Typically 
it is desirable to decompose the magnetic field B as 
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£ « <£> + 6 £ 


( 1 ) 


where <...> denotes an averaging operation and 6£ is the 
fluctuating field. This decomposition is of essential impor- 
tance in MHD because the mean field cannot be eliminated by 
transformation to a moving frame CElsasser, 19563 and it 
produces dynamical couplings that produce anisotropies in both 
wave-like C Cowling, 19763) and turbulent CShebalin, Matthaeus 
and Montgomery, 1983; Montgomery, 19833 magnetohydrodynamics. 

In the interplanetary medium, the mean field is not a DC field - 
there are time variations in the field at essentially all 
scales. On the other hand, MHD is local and allows small scale 
fluctuations to respond to the slowly varying long wavelength 
component of the magnetic field as though the magnetic field 
were strictly uniform and constant. In some cases, interplane- 
tary Intervals have unambiguous r.ean field values (e. g. , 
Alfvftnic periods of sufficient duration). In general, however, 
the decomposition into mean and fluctuating fields may be 
modified in several ways. First, the calculated value of the 
mean may not be stable. For example, the mean might undergo 
unacceptably large changes if the endpoints of the interval are 
shifted across a sector boundary. In addition, if time domain 
windowing or detrending is performed, the calculated mean would 
be changed. Irrespective of the implications of the dynamical 
interpretation, the variances and correlation times will bt ill- 
defined whenever the decomposition represented by equation (1) 
is ambiguous. 

The basis for alleviating some of the ambiguities in data 
selection can be found in the theory of stationary random 
processes CPugachiev, 1962; Panchev, 19713. A method has been 
developed to quantitatively test the hypothesis that interplan- 
etary magnetic field data intervals that are either of suffici- 
ent duration or otherwise properly chosen behave as if they are 
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weakly stationary random functions of time CMatthaeus and 
Goldstein, 1982b3. By definition, a weakly stationary magnetic 
field signal is one for which the ensemble values of the mean 
and che correlation functions are Independent of the selection 
of the origin of time. Under rather mild conditions, an ergodic 
theorem guarantees that time averages converge to ensemble 
averages as the duration of the data sample is increased without 
bound. It can be shown that this implies not only that the 
value of <B> calculated from the data converges to its ensemble 
value, but also that the form of the convergence is specified. 

In cases where convergence to weakly stationary behavior can 
be demonstrated, the mean field, variances, correlation func- 
tions and spectra can be reliably obtained. Matthaeus and 
Goldstein C1982b3 showed that the theoretical predictions for 
the rate of convergence of means and variances of magnetic field 
components agree remarkably well with the appropriate quantities 
obtained from solar wind data comprising many solar rotations. 
There are, however, several difficulties with simply using such 
very long data intervals to obtain highly convergent statistical 
estimates of ensemble properties of the solar wind. Long and 
continuous data sets are not always available. More fundamen- 
tally, the mean field which can he accurately obtained by 
averaging over many solar rotations is not the local mean field 
to which MHD fluctuations and plasma particles respond. In 
fact, averaging over many months causes a cancellation of the 
sector structure through many reversals of sign, and <B> tends 
to lie in the direction normal to the ecliptic plane. Statisti- 
cal quantities so obtained have lost information about the 
dynamically important local mean field and are biased by 
inappropriately classifying very low frequency changes of the 
solar field as fluctuations. Furthermore, long data Intervals 
disallow CMattheaus and Goldstein, 1982a; Goldstein, Burlaga and 
Matthaeus, 19843 the conversion from frequency to wavelength in 
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the super-Alfvdnlc solar wind flows. This is undesirable since 
turbulence theory usually emphasizes the role or spatial 
correlations and wavenumber spectra. 

What is needed are methods to calculate convergent statis- 
tics using intervals that are not so long as to wash out 
Information about the local mean field. Fortunately, an 
adequate degree of convergence to the predictions of ergodlc 
theory has also been obtained for Intervals which are not many 
solar rotations in duration, but are still long compared with 
the correlation time CMatthaeus and Goldstein, 19B2bD. In some 
cases, it Is possible to select data Intervals that do not 
contain individual anomalies such as sector crossings or shocks 
and subsequently verify that the theoretically expected conver- 
gence properties are well satisfied. Conversely, a "stationari- 
ty test” procedure can be constructed in which an arbitrary data 
set is compared to the predictions of the ergodic theory and the 
data is then judged as "stationary" or "nonstationary" in 
accordance with a quantitative figure of merit of agreement with 
the theory. 

In the following section we describe in detail a method for 
obtaining convergent estimates of ensemble properties of the 
interplanetary magnetic field. The basic assumption is that the 
collection of a large number of statistically stable estimates 
of the field properties constitutes an adequate representation 
of the ensemble. The principal data selection criterion will be 
the stationarity test method applied to data intervals for which 
a local mean field is physically meaningful and which do not 
violate the assumption of local spatial homogeneity. In this 
way we seek to characterize the interplanetary magnetic field in 
a manner that is relevant to statistical theories of MHD 
turbulence. 
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M ethod 


He assume that one or more sets of magnetic field data are 
available , which we designate as S* where 


S a = { B j ( t a + nAt) |n* 0,1... N a ; j - 1,2,3) (2) 

Indicating that three orthogonal components (j * 1,2,3) arc 
known for the magnetic field, g, at time values of t + nAt. 

Data, set S contains values of g. These are taken for 
convenience to have a constant time spacing At. The division of 
available data into sets labelled by different values of a Is 
somewhat arbitrary and has significance mainly as a computa- 
tional convenience. For example, a large data set with a data 
gap of significant duration in the middle should probably be 
broken into two sets. On the other hand, small data gaps can be 
included in a single set by the usual method of flagging "bad 
points" or data gaps. It is int'jnc'ed that each set S a be large 
enough that it would show highly convergent results in a 
statlonarlty test. However, if an interplanetary data set Is as 
long or longer than a solar rotation, then the test used on the 
entire set must use the form of the ergodic theory that includes 
solar rotation effects (see equation 12 in Matthaeus and 
Goldstein C1982b3. 

He do not perform statlonarlty tests on the entirety of each 
interval S a in view of the very long computational time 
required. However, a simple criterion for the acceptability cf 
a set S a is that it should span a very large number of 
correlation times, say, several hundred. Six months to & year 
of data is a convenient interval. It would be possible to 
proceed with sets that are not of such large duration, but then 
a larger number of sets would be required. In the analyses 


- 11 - 


described below, date gaps have been filled by linear 
Interpolation and a low band-pass digital filter was used to 
suppress aliasing. 

The main part of the technique consists of operations on 
each of the sets S a . An automated procedure is established. In 
tf ' fora of a modularized FORTRAN code to search the data in 
each set for subintervals of a predetermined length, and to 
store statistical information from those subintervals judged 
adequately stationary for later use. In advance of a single 
"run", parameters must be specified governing the selection of 
subintervals and determining the statlonarlty test criteria. 

The steps taken In each run, which are also depicted in Figure 
1, are: 


I. Select a subinterval from S a . For a given run, the 
number of points (time values) is fi^ed at NLEN, which should be. 
large enough to correspond to at least several correlation 
times. The next subinterval begins at a fixed number of points 
(N5LIDE) from the last interval's starting time index. The 
position within the set S a is analyzed to check if further 
processing Is requested. 

II. Preliminary statistics are calculated from the subin- 
terval dAta. First, the subinterval is rejected if the number 
of bad points flagged exceeds a predetermined percentage. 

Plasma velocity data is also used to calculate mean radial flow 
speeds in the subinterval for the purpose of converting frequen- 
cy data to wavenumber data. If the speed data has too many 
gaps, the mean solar wind speed defaults to a previously 
specified value (usually 400 kp/s), otherwise the mean speed is 
calculated. Next, the mean values and variances of the magnetic 
field components are calculated. 
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III. The data are transformed to a mean iieJd coordinate 
system. The elements of the rotation matrix are retained, as 
are the calculated means and variances of the components in the 

new coordinates. The data we used are initially in heliocentric 

/\ 

coordinates. In this coordinate system R is along the sun- 

A 

spacecraft line, positive away from the sun, T is perpendicular 

A AAA 

to R, and N = R x T. The rotated coordinates are in the 

A __ A 

directions B x N, B x (B x N) and B, where B is the mean field 
vector calculated from the subinterval. 

IV. A Blackman-Tukey mean-lagged-product correlation 
analysis is performed on the subinterval data using the method 
described by Matthaeus and Goldstein C1982a3. Estimates of the 
three diagonal and three antisymmetric off-diagonal elements of 
the correlation matrix CMatthaeus and Smith, 19813 are calcu- 
lated and retained, as are the statistical weights of each 
estimate and the values of the corresponding spatial lag. The 
latter is formed by multiplication of the time lags by the mean 
solar wind speed in the subinterval. The correlation time and 
length are also calculated. If the subinterval contains less 
than NC correlation times, the subinterval is rejected (usually 
we take NC = 10). Note that even though the mean field direc- 
tion varies according to the particular subinterval, the spatial 
lags always correspond to radial separations. Information about 
the statistical weight of each correlation estimate can be 
easily used later to calculate composite correlation functions. 

V. A stationarity test is performed on the three components 
of the magnetic field in the subinterval. The procedure 
consists of averaging the field data over time intervals of 
duration x, where x varies from the sampling time up to one- 
fifth of the subinterval duration, NLEN At. We require that the 
variance of these averages calculated from all such samples 
within the subinterval be sufficiently close to the value 
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expected for a stationary random process converging to its 
ensemble mean. If an average over time x of field component 
is denoted by CBj3 T and the variance of such an average within 
the subinterval is denoted by A 2 CBj3 T , then the relevant 
convergence behavior is given by CMatthaeus and Goldstein, 

1982b, equation 111 

A 2 CB .3 = 2(x ) . a. 2 It (3) 

3 t c 3 3 

where (t ) . and a 2 are the correlation time and variance of the 
c 3 J 

jth component calculated for the entire subinterval in step IV. 
The weighted vector version of (3), 


A 2 CB3 la 2 = > A 2 CB .3 / > a 2 (4) 

~ £_j J ^ il_| J 

j j 

is used in the stationarity test Ccf. Matthaeus and Goldstein, 
1982b3. 

The correlation analysis is done first so that these 

parameters are available for the stationarity test. The left 

hand side of (4) is calculated from the subinterval data and 

compared with the right hand side for specified values of t c /t. 

We have usually chosen x /x values of one and two. 

c 

The subinterval is judged as adequately stationary if the 

fractional difference of the empirical value of the left hand 

side of (3) and the right hand side is less than one at x c /x = 1 

and if it is also less than one-half at x /x = 2. These 

c 

parameters control the sensitivity of the test and are specified 
at the beginning of the run. If the subinterval is judged 
"nonstationary", it is rejected. 
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VI. The statistics and characterizing parameters of 
stationary intervals are stored. This includes means, 
variances, correlation lengths and correlation functions, as 
well as the position of the subinterval within S a , the mean 
solar wind speed, the rotation matrix, the correlation lags and 
the statistical weights of the correlation estimates. Process- 
ing then continues with step 1. 

Steps I through VI are applied to all selected sets S ot . 
Several such runs through all the data can be done with differ- 
ing values of NLEN, and with various values of stationarity test 
sensitivity. The entire procedure (see Fig. 1) results in a 
collection of statistical analyses that may be quite large in 
number, all of which have been calculated from a method involv- 
ing minimal subjectivity. The data base obtained may then be 
treated as a representation of the magnetic field ensemble 
corresponding to the conditions prevalent in the data intervals 
S* 


ISEE data: A Magnetic Ensemble at 1 AU . 

The method has been used to analyze ISEE-3 magnetic field 
data made available to us under the auspices of the ISEE Guest 
Investigator program. We are indebted to E. J. Smith and B. 
Tsurutani for providing us with the data that made this study 
possible. At various stages we have made use of ISEE plasma 
velocity data (P. I.: S. J. Bame, LAND from the National Space 
Science Data Center (NSSDC) Omnitape and ISEE-3 helium speed 
data (provided by K. W. Ogilvie) to estimate mean solar wind 
speeds. 

We have used two time intervals of magnetic field data. The 
five minute data samples were linearly interpolated over gaps, 
filtered with a 128 point Remez finite impulse response filter 
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[McClellan, Parks and Rabiner, 19793, and decimated to every 
third point. This preprocessing produced fifteen minute 
averages relatively free of aliasing errors. The first data 
set, S 1 , consists of 27,736 fifteen minute averages, spanning 
the interval from day 152 (1978) to day 228 (1979). The second 
set, S 2 , begins on day 180 (1979) and ends on day 365 (1979), 
for a total of 17,814 samples of each component. The total 
percentage of gaps and bad points was 11.5% in S 1 and 11.9% in 

n2 


Runs were performed on S 1 and S 2 according to the procedure 
given in the last section. Twelve runs were made using six 
subinterval lengths on each set. We selected subinterval 
lengths (NLEN in step I above) of 60, 125, 250, 500, 1000 and 
1500 fifteen minute points. The criteria for acceptability of a 
subinterval were, for all cases: (i) at least ten correlation 
times are present in the subinterval and (ii) the stationarity 
test was satisfied in that the data was consistent with eq. (3) 
evaluated at the points and with the precision described in step 
V above. Criterion (i) must be satisfied to perform the 
stationarity test. The parameter for incrementing the starting 
point when selecting a new subinterval (NSLIDE in step I) was 
either 50 or 100 in all twelve runs. 

Table 1 summarizes the run parameters, and indicates the 
number of acceptably stationary intervals found in each case. 

The total number of data intervals considered was approximately 
355 for each of the six set S 2 runs, and was approximately 550 
for the set S 1 runs with NSLIDE = 50, and half that number for 
the set S 1 , NSLIDE = 100, runs. Using this information, a 
percentage of intervals found to be stationary is included in 
Table 1. The percentages found for both S 1 and S 2 when NLEN = 
250 and 500 are quite similar. These subinterval lengths 
correspond to durations of 2.6 and 5.2 days. The frequency of 
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finding stationary data subintervals falls off by about 10% when 
NLEN is increased to 1000, corresponding to subintervals of 10.4 
days duration. The NLEN = 1500 runs (15.6 days) are found to be 
stationary even more rarely. 

The relative probability that data sets of a few days 
duration are stationary is consistent with the viewpoint 
CMatthaeus and Goldstein, 1982a3 that spectral analysis of 
stationary data requires samples of at least several correlation 
times. At first sight, the decrease in likelihood that longer 
data sets are stationary might appear puzzling because one 
expects that longer samples should behave mere like the ensem- 
ble, and therefore should be more consistent with eq. (3). 
However, it must be borne in mind that the ensemble we are 
constructing is cne that is consistent with the notion of 
stationarity about a local, well defined mean magnetic field 
direction. We are separating the behavior with respect to the 
local mean field from the behavior due to the systematic 
variations in the field produced by phenomena linked to the 
solar rotation period (see discussion in the paragraph following 
equation (1)). The use of eq. (3) above rather than eq. (12) in 
Matthaeus and Goldstein C1982b3 (which includes solar rotation 
effects) precludes finer agreement in the longer data subinter- 
vals as effects of solar rotation become significant. The NLEN 
= 1500 runs are more than half a solar rotation period, and the 
decrease in percentage of stationary intervals for those cases 
simply reflects the dependence of the mean field calculated over 
that interval on the systematic behavior of low frequency 
structures in the solar wind. 

We now turn to a discussion of the some of the statistical 
information that can be gleaned from the output of these runs. 

Average Properties of the Ensemble . 
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Averages may be calculated from the data produced by the 
runs listed in Table 1 in several meaningful ways. A statisti- 
cal quantity calculated in each of the stationary subintervals 
may be averaged over ail of the 1,959 data sets produced without 
regard for subinterval length. Averages may also be accumulated 
separately from each of the twelve runs listed in Table 1. 
Possible dependence of the averages on the subinterval length 
can be investigated by forming averages using all subintervals 
having a specific length in both S 1 and S z . Finally, just che 
data in S 1 or in S 2 may be used, again ignoring subinterval 
length. We will associate the symbol "S tot " with averages 
accumulated from all the stationary subintervals in the union of 
S 1 and S z . When an average is calculated from all of the 60, 
125, 250, 500, 1000 or 1500 point stationary subintervals in the 
union of S 1 and S 2 , it will be associated with the symbol " s 60 ” » 

" S z,o"' " S »oo"' " s ,ooo"' or " s i,oo" respectively. To 
avoid confusion, an overbar will indicate averaging within a 

subinterval, while angle brackets will indicate averaging over a 
collection of subintervals. 

The first class of averages we discuss pertains to mean 
field values in the subintervals. Table 2 shows averages of the 
means of the R, T and N components of the magnetic field, the 
average vector magnitude, the average of the angle 0 between the 
radial direction and the mean field, and the average deviation 
of 0 from 90°. These quantities are given for the entire 
collection of stationary subintervals, S tQt , for the 60 point 
subintervals, and for the 1500 point subintervals, S 

The first thing to notice is that the average magnitude is 
larger for the R and T components than for the N component, in 
qualitative agreement with the Parker model of the mean inter- 
planetary field. The average of the mean vector magnitude is 
5.3 nT when calculated from S go . This is comparable to quoted 
values of typical low speed "structureless" solar wind CHund- 
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hausen, 1972, page 443. However, values of the average magni- 
tude of the local mean calculated from both S. . and S,_ are 
somewhat smaller, which probably arises because low frequency 
fluctuations are included. He have not included the contribu- 
tion to the average field magnitude due to the fluctuations in 
the subintervals, i. e. we have calculated the magnitude of the 
average field, not the average of the magnitude of the field. 

The average of 0 is 97° indicating nearly complete cancella- 
tion of inward and outward directed fields due to the sector 
structure. The average deviation of 0 from 90° is about 38° for 
S tot' an< * systematically decreases from 40° in S gQ down to 30° 
for S averages. This suggests that the average mean fields 

1900 

subtend an angle of about 52° with respect to the radial 
direction. Note that this angle includes contributions from the 
small but non-negligible normal component of the mean field. 

The mean solar wind speed is 405 km/s for the entire ensemble 
and ranges from 411 km/s in down to 394 km/s in the 
averages. This indicates a slight systematic correlation of 
higher solar wind speed with stationary periods of shorter 
duration. High speed flows generally tend to be shorter than 
15.6 days (the duration of a 1500 point subinterval). Similar- 
ly, these results suggest that low speed flows are more likely 
to be stationary over longer time periods than high speed flows. 

The mean solar wind speed in set S 1 was 377 km/s compared 
with 429 km/s for set S 2 . This effect is larger than the 
variation of speed with subinterval length and may be a conse- 
quence of the fact that set S 2 covers a time interval closer to 
solar maximum than set S 1 . However, this conclusion is uncer- 
tain in view of the relatively high percentage of data gaps in 
the solar wind speed data used. 
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Statistics involving the variances of the magnetic field are 
presented in Table 3. Ensemble averages of the variances of the 
field components in heliocentric and in mean field coordinates 
are given, as are the ensemble average values of total field 
variance, El 2 , the fluctuation magnitude, SB = J ( 6B- SB) , and the 

value of <8B/|B|>. These data are shown for S. . , S„ . S , 

S_ and S . Recall that the ensemble averaging procedure 
does not commute with nonlinear algebraic operations on field 
quantities in each subinterval. Thus, for example, <<SB/|5|> is 
not necessarily equal to <SB>/<|B|>. 

One interesting feature of this data is that the field 
variance is quite large for the component of the fluctuations 
parallel to the mean field. The ratio is about 8:9:10 for the 
entire ensemble and is about 11:16:10 for the subintervals. 
The dominance of the fluctuations parallel to the mean field 
increases up to a ratio of about 6:7:10 for the long S 15Q0 
samples . 

The total field variance systematically increases as the 

subinterval size is increased from about 25 nT 2 for S samples 

60 

to 55 nT 2 for S, samples. The quantity <SB/|B|> is about 2.3 
for the entire ensemble, but also shows systematic dependence on 
the sample size, varying from 0.95 for the shortest subintervals 
up to 5.0 for the longest. This behavior of both the variances 
and the values of <SB/|B|> supports the view that long- 
wavelength fluctuations play an important role in determining 
the structure of the interplanetary field CBurlaga, 1983; 
Goldstein, Burlaga and Matthaeus, 19843 and may also have 
implications for theoretical treatment of particle motions in 
the heliosphere. This issue will be taken up in the discussion 
section. 
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Averages of the correlation lengths and correlation times of 
the components and the total field are shown in Table 4 for 
S. . , S a . and s , „ „ . samples. The value of the 

correlation time is calculated from the correlation function 
CMatthaeus and Goldstein, 1982a3 for each subinterval. Frozen- 
in-flow CTaylor, 19383 is used to convert correlation times to 
correlation lengths usi >g the mean solar wind speed in the 
subinterval. Again there, is a systematic increase in the 
correlation lengths and tin. •*s as the subinterval size increases. 
The total ensemble value of the correlation length L c is 4.9 x 
10 11 cm; the value is 1.2 x 10 11 cm and the S value is 
1.5 x 10 1 cm. There is also an anisotropy in the correlation 
lengths in all but the S gQ run with the largest values appearing 
for the components parallel to the mean field. The correlation 
length in the S go run is almost perfectly isotropic. The same 
pattern of anisotropy and dependence on interval size is seen in 
the correlation time averages, which indicates that there is not 
a strong positive correlation between increasing solar wind 
speed and correlation time. The correlation time is about 3.4 
hours for S tQt , and ranges from 0.81 hours for S g0 up to 10.4 
hours for S _ . In all cases the average correlation time is 

1500 

considerably less than the duration of the subintervals (see 
Table 1) which is a posteriori a requirement for accurate 
determination of T c . (The expression we use for T c guarantees 
that T c becomes comparable to NLEN At if all the power in the 
signal is at the lowest sampled frequency Csee Matthaeus and 
Goldstein, 1982a3.) 

Other systematic variations in the statistics can be 
investigated by calculating separately averages from the S 1 and 
S z subintervals. Since S 1 and S 2 are data sets of comparable 
length, comparisons of the average s-atistics from these 
separate analyses provide a certain measure of confidence that 
the ensemble average results represent a fairly stable statisti- 
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cal representation of the solar wind, at least during this 
portion of the solar cycle. 

In Table 5 the correlation lengths and times and variances 
of the components and the total field, and the magnitude oi the 
local mean field components are compared for the major intervals 
S 1 and S 2 . The entries in Table 5 should be compared with the 
corresponding entries in Tables 2, 3, and 4. The values of the 
correlation length calculated from set S 1 differ by just a few 
percent from those calculated from set S z . For example, the 
total correlation length from set S 1 , 5.12 x 10 12 cm, differs 
from the set S 2 value, 4.73 x io 12 cm, by 8%. The correlation 
times are actually more variable than the correlation lengths, 
with differences in S 1 and S 2 estimates of the total correlation 
time differing by 20%. The variances of the field calculated 
from S 1 and S 2 also compare well. The total variance from set 
S 1 is 40.0 nT 2 and from set S 2 is 41.0 nT 2 which differ by only 
2.5%. In contrast, the average magnitude of the local mean 
field varied by 14% in the two cases (3.77 nT in S 1 and 4.39 nT 
in S 2 ). Evidently the statistical averages are reproducible 
with good precision, particularly the correlation lengths and 
field variances which are closely tied to the structure of the 
turbulence. Nevertheless, the standard deviations of these same 
statistical estimates, shown in Tables 2, 3, and 4, indicate 
that there is a considerable spread in the estimates in the 
subintervals. For example, the standard deviation of the 
estimates of the total field variance from the entire set S fcot 
was 34 2 nT 2 (Table 3), a significant fraction of the average 
variance, 40.6 nT z , even though the average variance was almost 
identical in S 1 and S 2 . 
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Distribution of Values in the Ensemble. 


Each of the statistical quantities in Tables 2-5 may be 
viewed as properties of a random variable with well-defined mean 
value. The distribution of these variables in the ensemble can 
be studied by calculating the frequency of occurrence of values, 
which is equivalent to a coarse-grained zero-lag probability 
distribution function. Extensive correlation analysis would 
contain information characterizing more complete probability 
descriptions of these random variables. A complete presentation 
of these properties will not be undertaken here. Rather, we 
shall focus on a small subset of diagnostics that describe 
distributions which may be of broad interest. 

Figure 2 shows a histogram of solar wind speed values for 
the entire ensemble. In this and all subsequent figures, 1,958 
data sets are used; i. e., one data set with an anomalously high 
variance (in excess of 525 nT 2 ) was dropped from the ensemble. 
The distribution is peaked just below the mean speed of 405 
km/s, and has a high-speed tail extending beyond 600 km/s. The 
shape of the distribution and its mean value are quite 
comparable to the results obtained from Vela 3 during 1965 to 
1967 CHundhausen, 1972, page 44; Hundhausen et al., 19703. 

A histogram of values of the S fcot magnetic field variance, 
SB 2 , is given in Figure 3. The distribution again displays a 

tail at large values of the variance, so that the mean value 

lies a little above the main peak of the distribution. There is 

also a rapid fall-off just above 60 nT 2 . 

Figure 4 is a histogram of the correlation length calculated 
from the S fcot data. The main maximum of the distribution is 
again at a value somewhat lower than the mean of 4.9 x 10 11 cm. 
The maximum is displaced from the mean by a relatively greater 
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amount than In Figure 3 because the h\gh correlation length tall 
Is more populated than the high variance tall. In fact, there 
Is a small but noticeable secondary peak In the population near 
2.3 x 10 12 cm which Is predominantly due to the 3 and S 

1000 1900 

data Intervals. 

Scatter plots are convenient for displaying qualitative 
Information about she elements of the ensemble. For each 
stationary subinterval a point is plotted in a plane whose axes 
are each a characterizing statistical parameter. Figure 5 Is a 
scatter plot with the variance plotted on the ordinate and the 
solar wind speed plotted on the abscissa. A total of 1,958 
points are plotted. Hhlle there are a few points that have 
large variance in conjunction with high speed, the overall 
impression is that there ib little correlation between the 
variance and mean solar wind speed in the collection of inter- 
vals Sj. Q £. The small number of samples at speeds greater than 
575 km/s prohibits drawing firm conclusions about possible 
correlation between very high speed and very high variance. 

A scatter plot of the total correlation length vs. solar 

wind speed is shown in Figure 6. The great majority of points 

lie below 10 12 cm, and most of these are for intervals in which 

the speed is less tnan 400 to 450 km/s, which also follows by 

inspection of Figures 3 and 4. There is a subpopulation of 

intervals with L on the order of 10 12 cm and solar wind speed 

c 

less than 400 km/s. There is an even more sparse subpopulation 
at L c ranging from 1.5 x 10 12 cm up to 3 x 10 12 cm. These seem 
to be less concentrated at speeds under 400 km/s. The latter 
subpopulation is associated predominantly with the longer data 
intervals S iooQ and S 150Q . The shorter intervals, such as S 60 , 
have lower L c and slightly higher than average speed. 
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Finally, Figure 7 is a scatter plot of the correlation 

length vs. the field variance. Here we see a concentration of 

Intervals with variance below 50 nT 2 and L near a few tines 

c 

10 11 cm. This region is preferentially populated by the shorter 
intervals and also appears to have a broad but discernible 
positive correlation between variance and correlation length. 

In the area around 5! 2 = 50 nT 2 the points span a wide range of 
values of L c , but are still mostly below 10 12 cm. Higher 
variance intervals are also mostly found with L c < 10 12 cm, but 
several patches of points are found at higher L c> 

Discussion . 

The method we have used to characterize the statistical 
properties of the interplanetary magnetic field involves 
averaging over many subintervals of available data without 
regard to classification into specific types of intervals such 
as "transient flows" or "corotating flows" Ccf. Goldstein, 
Burlaga, and Matthaeus, 19841. The features of the data we have 
emphasized are the local "stationarity property", stability of 
calculated local mean magnetic field, and effects of subinterval 
length. 

One point of contrast between previous results and those 
presented here concerns the issue of "Alfv6nic" periods. Since 
the work of Coleman C1968J, Belcher and Davis C19713 and others, 
a considerable amount of effort has been directed towards 
investigating interplanetary data intervals that show a high 
degree of correlation between fluctuating magnetic and plasma 
velocity fields. Belcher and Davis found that these periods are 
characterized by fluctuations whose variance is least in 
component parallel to the local mean magnetic field, 1. In the 
wave picture of MHD, this corresponds to a spectrum of 
fluctuations with wavevectors predominantly parallel to 1, i. e. 
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"slab" geometry, and has led to use of "slab" models for the 
scattering of cosmic rays In the Interplanetary medium CFisk, 
19793. A comprehensive survey of the frequency of appearance of 
Alfvdnlc periods would be necessary to judge the extent to which 
characterizing interplanetary fluctuations as being Alfvdnlc is 
appropriate. 

Khile the present study has not attempted to Isolate 
Alfvdnlc periods, our results do bear on this Issue. The value 
of the ratio of variances In Alfvftnic periods is found by 
Belcner and Davis C19733 is 5:4:1. In contrast, we find that 
the ratio Is 8:9:10 for the total ensemble, 6:7:10 for the 1500 
poin.. data intervals, and 11:16:10 for the 60 point data 
Intervals. In all subinterval sets longer than and including 
the 500 point set, the maximum variance was found to be in the 

component parallel to j£. In the S , S and sets, the 

maximum variance was in a component other than the parallel 
component, but only in the S set was the minimum variance 

6 O 

found In the parallel component. Even in the S run, the 

6 O 

variance ratio was quite different from that found by Belcher 
and Davis C19713 to be typical of Alfvenic periods. This 
strongly suggests that it is inappropriate to associate the 
properties of Alfvtnlc periods with interplanetary turbulence in 
general. For all interval sizes used, the results obtained here 
give component variances which are much closer to Isotropic than 
slab- like. 

The calculated average correlation lengths of the components 
(see Table 4) are In the approximate ratio 8:6:10 for the entire 
ensemb.lc , but they also show a systematic dependence on Interval 
length. In particular, the degree of anisotropy of component 
correlation lengths decreases with decreasing data Interval 
size. The correlation lengths of the three components are 
essentially identical fcr the S 6Q data intervals (Table 4). 
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This does not imply that the wavevectors are isotropically 
distributed, but it does Indicate that certain anisotropic 
features are more pronounced in the very long wavelength part of 
the fluctuation spectrum. 

One of the most significant features of these results is 
that the variances, correlation lengths and values of SE/|B| 
systematically depend on the data interval length (see Tables 3 
and 4). As the interval length is increased, lower frequency 
signals are sampled, presumably including signals associated 
with stream and sector structure and other large scale coherent 
features of the heliosphere. There is also the possibility that 
the random component of the turbulence contains "excess" energy 
at very low wavenumber, due to self organization effects such as 
the inverse cascade of magnetic helicity CFrisch et al., 19753 
or selective decay processes [Montgomery, Turner and Vahala, 
1978; Natthaeus and Montgomery, 19803. Some indication of this 
low wavenumber behavior has been seen in analysis of Voyager 
data [Goldstein, Burlaga and Matthaeus, 1984; Burlaga and 
Goldstein, 19843. 

The present results support the view that the very large 
scale structures in the interplanetary medium affect statistical 
estimates in the following way: The variance of the magnetic 

field, which is proportional to the local mean magnetic energy 
density of the fluctuations, increases monotonically as the data 
interval size is increased from 0.625 days to more than 15 days. 
Accordingly, estimates of the strength of the mean magnetic 
field decrease so that the calculated values of Sl/|B| decrease 
as the interval size increases. Similarly, che correlation 
lengths and correlation times of the components and the total 
field also monotonically increase as the interval size is 
increased over the range we have examined. 
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Estimates of the correlation time calculated from intervals 
much longer than those used here CMatthaeus and Goldstein, 

1982b; Goldstein, Burlaga and Matthaeus, 19843 are only slightly 
greater than the value we find for the longest intervals. For 
example 15 hours is a typical correlation time calculated from 
intervals ranging from two to 22 solar rotations while the S 

1500 

estimate above corresponding to about a half solar rotation is 
about 10 hours. This indicates that the changes seen above as 
the interval size is increased may saturate at times comparable 
to the solar rotation period. However, the method used here 
becomes increasingly ineffective for interval lengths longer 
than 15 days because the frequency of finding stationary 
intervals decreases rapidly for longer intervals (see Table 1). 
We attempted a 3000 point run on set S 1 and found no stationary 
intervals. 


Summary and Future Directions. 

In this paper we have presented a technique for extracting 
high statistical weight statistical estimates from a collection 
of data intervals. We have applied the method to the study of 
interplanetary magnetic fluctuations at 1 AU using ISEE-3 data. 
The focus of this method is the characterization of interplane- 
tary fluctuations in the context of a locally homogeneous, time- 
stationary, MHD turbulence description. The local mean magnetic 
field plays a central role in such a theoretical model. 

Therefore the method has been constructed to choose intervals 
that are judged to be weakly stationary with respect to a local 
mean field in accordance with ergodic theory. No further data 
selection criteria were used. The output of the computational 
algorithm consists of a large number of data sets each of which 
contains statistical information from a single stationary 
subinterval. Further analysis was presented giving the ensemble 
average values of mean and fluctuating field statistics and 
correlation lengths and times. 
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The main results found from this preliminary study indicate 
that the correlation times and lengths and the field variances 
and estimates of mean field strength have reproducible values 
that nonetheless have a systematic dependence on the length of 
the data subinterval. We also found that the variances of the 
components of the magnetic field calculated from the ensemble do 
not have the characteristics of Alfv6nic periods. 

The method presented may find future application in the 
study of data from other spacecraft when long intervals of 
nearly continuous data coverage are available. Similar tech- 
niques can be applied to plasma velocities , densities , and other 
experimental data in turbulent plasmas. Additional physical 
problems of the interplanetary magnetic field at 1 AU can also 
be studied with this method. For example, the one-dimensional 
magnetic correlation functions can be calculated from the 
ensemble. Furthermore, because the ensemble has retained 
information about the direction of the mean magnetic field in 
each subinterval, it is possible to calculate a three-dimension- 
al correlation function and spectrum from the ensemble under the 
assumption that the turbulence is axisymmetric with respect to 
the local mean field direction. It is likely that the ensemble 
methods can be generalized to calculate the one-point probabili- 
ty distribution function of the magnetic field fluctuations, 
extending the work of Whang C19773. Finally, in a separate 
publication, we will report on a study of the dependence of the 
variance and correlation times on the data interval length which 
may further elucidate the role of large scale structures in the 
solar wind. 
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TABLE 1. Stationary Subintervals In Sets S 1 and S z 


NLEN 

(15m 

NSLIDE 
samples ) 

DURATION 
(days ) 

NUMBER 

STATIONARY 

PERCENT 

STATIONARY 

(set S 1 ) 

60 

50 

0.625 

231 

42% 

125 

50 

1.3 

251 

46% 

250 

50 

2.6 

250 

45% 

500 

100 

5.2 

119 

43% 

1000 

100 

10.4 

94 

34% 

1500 

50 

15.6 

98 

18% 

(set S 2 ) 

60 

50 

0.625 

166 

46% 

125 

50 

1.3 

170 

48% 

250 

50 

2.6 

172 

48% 

500 

50 

5.2 

177 

50% 

1000 

50 

10.4 

132 

37% 

1500 

50 

15.6 

99 

28% 
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TABLE 2. Ensemble Averages of the Mean Field 
(and Standard Deviations) 




S 

tot 

S 

60 

S 

1900 

Hr' 

(nT) 

2.44 

(1.6) 

3.21 

(1.8) 

1.06 

(0.8) 

|BJl 

(nT) 

2.80 

(1.7) 

3.49 

(2.0) 

1.68 

(1.0) 

* ®N * 

(nT) 

0.90 

(0.9) 

1.31 

(1.3) 

0.42 

(0.3) 

111 

(nT) 

4.11 

(2.1) 

5.30 

(2.2) 

2.16 

(1.1) 

0 

( °) 

96.6 

(41. ) 

96.7 

(43. ) 

99.8 

(34. ) 

|0-9O° | 

(°) 

35.9 


38.9 


30.4 


V 

( km / s ) 

401. 

(69. ) 

411. 

(79. ) 

394. 

(47. ) 

sw 
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TABLE 3. Ensemble Average Variances of the Magnetic Field 

(and Standard Deviations) 



S T0T 


o 

10 

CO 

S 

2 S 0 

S 500 

S 1 500 



11.4 

(12) 

5.7 

(8) 

10.9 

(10) 

13.0 

(8) 

15.2 

(7) 

si; 2 

16.5 

(16) 

9.0 

( 18) 

16.2 

( 16) 

21.3 

(12) 

25.2 

(9) 


13.7 

(12) 

10.3 

(15) 

14.5 

(13) 

16.3 

( 10) 

15.1 

(7) 

sir 2 

11.4 

(12) 

7.5 

( 13) 

12.3 

(14) 

14.0 

(11) 

14.3 

(9) 

si; 2 

14.1 

(13) 

10.7 

(17) 

14.9 

(12) 

16.7 

(9) 

16.5 

(8) 

si; 2 

15.0 

(15) 

6.9 

(14) 

14.4 

(15) 

19.9 

(13) 

24.7 

( 10) 

5B 2 

40.6 

(34) 

25.0 

(39) 

41.6 

(35) 

50.6 

(27) 

55.5 

(18) 

51 

5.9 

(2) 

4.4 

(2) 

6.03(2) 

6.89(2) 

7.36 

(1) 

51/ 1 B I 

2.26 

(3) 

0 . 95 ( .7) 

1.69(1.4) 

2.84(2.6) 

5.02 

(4) 






TABLE 4. Ensemble Average Correlation Times and Lengths 

(and Standard Deviations) 





< T C 

In hours , 

L„ in 
c 

10 11 

cm) 




s 

tot 

s 


S 


S 


S 





60 

2 5 0 

900 

1900 

T 

Cl 

3.2 

(4.0) 

0.80 

( .15) 

2.2 

(1.0) 

3.9 

(1.9) 

9.6 

(7.4) 

T 

C2 

2.2 

(2.3) 

0.75 

( .22) 

1.8 

( .84) 

3.0 

(1.7) 

5.0 

(4.4) 

n 

*C3 

3.8 

(4.7) 

0.76 

( .24) 

2.4 

(1.0) 

4.6 

(2.0) 

11.8 

(8.2) 

T c 

3.4 

(3.7) 

0.81 

( .15) 

2.2 

( .68) 

4.1 

(1.3) 

10.4 

(5.7) 

L d 

4.6 

(5.6) 

1.2 

( .39) 

3.3 

(1.7) 

5.7 

(3.0) 

13.5 

(10. ) 


3.2 

(3.6) 

1.1 

( .39) 

2.7 

(1.3) 

4.3 

(2.5) 

7.2 

(7.3) 

L C3 

5.4 

(6.9) 

1.2 

( .31) 

3.5 

(1.7) 

6.6 

(2.9) 

16.8 

(12.1) 

K 

4.9 

(5.4) 

1.2 

( .31) 

3.3 

(1.2) 

5.9 

(2.0) 

14.8 

(8.5) 
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TABLE 5. Comparison of S 1 and S z 


S 2 


|BJl 

<nT) 

2.69 

2.69 

|B^| 

(nT) 

2.64 

2.95 

1 1 

(nT) 

0.87 

0.92 

"T 

(nT 2 ) 

9.5 

11.1 

ss: 2 

X 

(nT 2 ) 

17.2 

15.9 


( nT 2 ) 

13.3 

14.0 

SB 2 

(nT 2 ) 

40.0 

41.0 

T 

cx 

(hours ) 

3.6 

2.9 

T 

C2 

(hours ) 

2.5 

2.0 

T c, 

(hours ) 

4.3 

3.3 

T c 

(hours ) 

3.84 

3.06 

L c. 

( 10 1 1 cm) 

4.8 

4.6 

L c 2 

( 10 1 1 cm) 

3.4 

3.1 

L 

C 3 

( 10 1 x cm) 

5.7 

5.2 

L 

c 

( 10 1 1 cm) 

5.12 

4.73 
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Figure 1. Overview of the general algorithm uaed in 
constructing the Magnetic field enaeable. 
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Speed is plotted in Iw/s. 





Figure 3. A histogram of values of the S magnetic field variance 
(in nT 2 ). 
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Figure 5. A scatter plot with the variance (in nT ) plotted on 
ordinate and the solar wind speed (in fca/s) plotted on the abscissa, 
total of 1,958 points froa S fc t are plotted. 
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Figure 7. A scatter plot of the correlation length (divided by 10 
vs. the field variance (in nT 2 ). 
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